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ABSTRACT 

The orbital motions of halo stars in the Milky Way reflect the orbital motions of the progenitor systems in 
which they formed, making it possible to trace the mass-assembly history of the Galaxy. Direct measurement 
of three-dimensional velocities, based on accurate proper motions and line-of-sight velocities, has revealed that 
the majority of halo stars in the inner-halo region move on eccentric orbits. However, our understanding of the 
motions of distant, in-situ halo-star samples is still limited, due to the lack of accurate proper motions for these 
stars. Here we explore a model-independent analysis of the line-of-sight velocities and spatial distribution 
of a recent sample of 1865 carefully selected halo blue horizontal-branch (BHB) stars within 30 kpc of the 
Galactic center. We find that the mean rotational velocity of the very metal -poor ([Fe/H] < -2.0) BHB stars 
significantly lags behind that of the relatively more metal-rich ([Fe/H] > -2.0) BHB stars. We also find that 
the relatively more metal-rich BHB stars are dominated by stars with eccentric orbits, as previously observed 
for other stellar samples in the inner-halo region. By contrast, the very metal-poor BHB stars are dominated by 
stars on rounder, lower-eccentricity orbits. Our results indicate that the motion of the progenitor systems of the 
Milky Way that contributed to the stellar populations found within 30 kpc correlates directly with their metal 
abundance, which may be related to their physical properties such as gas fractions. The se results are c onsistent 
with the existence of an inner/outer halo structure for the halo system, as advocated by iCarollo et al.l 

Subject headings: Galaxy: evolution — Galaxy: formation — Galaxy: halo — Galaxy: kinematics and 
dynamics 



1. INTRODUCTION 

The luminosity of the Milky Way is dominated by its disk, 
where the great majority of stars (> 90%) are found. By com- 
parison, out of ~ 10 11 stars in the Milky Way, the stellar halo 
comprises only a tiny fraction (~ 1%), but this component is 
a precious source of information on the formation history of 
the Galaxy. First, halo stars are very old (~ 10-13 Gyrs), 
and their chemical compositions provide information on the 
ancient environments in which these stars formed. Typically, 
the metal ab undances of halo star s are less than 1/1 0th of th e 
Solar value (IChiba & Beersll2000t ICarollo et al.ll2007l 12010b . 
which immediately suggests that the metal enrichment due to 
supernova explosions in the Universe had not progressed very 
far when halo stars were formed. Secondly, the stellar halo is 
a collisionless system, hence two-body relaxation is expected 
to be unable to fully erase the initial orbital properties of halo 
stars. It follows that the present motions of halo stars reflects 
their motions in the early Universe, which can be used to ex- 
plore the kinematics of their progenitor systems, such as the 
gas clouds or dwarf galaxies in which these halo stars formed. 

To date, detailed analyses of the chemical and dynami- 
cal properties of halo stars has been confined to stars up to 
~ 10-15 kpc from the Sun (centered at the Galactocentric 
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distance of the Sun, ^8.5 kpc), although the spatial distribu- 
tion of the stellar halo extends to ~ 100 kpc or more. This is 
mainly because we do not possess sufficiently accurate mea- 
surements of proper motions for more distant halo stars. In 
order to avoid this limitation, many authors have studied the 
line-of-sight velocities and spatial distribution of halo stars, 
making use of various kinematic models - such as distribution 
function models of the stellar halo, or gravitational potential 
models of the Milky Way. However, the conclusions of pre- 
vious studies on the orbital distribution of halo stars well out- 
side the local region are divergen t. Some suggest tangentially- 
anisotropic orbita l distributions dSommer-Larsen et al.lll997t 
iKafle et ail 120 121), others su ggest radially-anisotropic distri- 
butions (IDeason et alj 120121). and still others suggest nearl y 
isotropic distributions (ISirko et all 120041 iThom et al.l 120051) . 
This might imply that the stellar halo is not a simple en- 
tity which can be described by a single distribution func- 
tion model. Indeed, b ased on observations of re lative nearby 
(d < 4 kpc) halo stars, (ICarollo et al.l2007l l2010b suggest that 
a dual-halo model is more appropriate, in which the stellar 
halo consists of a relatively metal-rich inner-halo component 
with a net zero to slightly prograde rotation, and a very metal- 
poor outer-halo component with a net retrograde rotation. Re- 
cent observations of retrograde outer-halo RR Lyrae stars by 
iKinman et al.l (120 12|). and a model -fitting analysis for distant 
halo stars by IDeason et al.l (1201 11) supports this idea, as do 
recent numerical s imulations of the form ation of Milky Way- 
like galaxies (e.g.. lMcCarthv et"aT1l2012l) . Here we introduce 
a new analysis of the halo system, which requires only a min- 
imum of assumptions, and does not require any kinematic 
models. 

This paper is outlined as follows. In Section 2 we describe 
our sample selection and expected errors in distance, line- 
of-sight velocities, and metallicities. Section 3 describes our 
analysis approach. In Section 4, we examine the application 
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of this technique to our sample of halo BHB stars, segregated 
on metallicity. Section 5 presents a brief discussion and con- 
clusions. 

2. SAMPLE SELECTION 

Our sample comprises 1865 blue horizontal-branch (BHB) 
stars from Data Relea se 8 of the Sloan Digital Sky Sur- 
vey (lAihara et al.ll201 ll) . with Galactocentric distances in the 
range 6 < r/kpc < 30, as or iginally selected and carefully val- 
idated by IXue et al.l d201 ll) . This sample is free from signif- 
icant contamination by the Sagittarius stream and the thick 
disk, as we apply a spatial masking scheme in their selection. 
Distances from the Sun are accurate to ~ 5- 10 %, and the 
line-o f-sight velocity errors are 5-20 km s -1 (see IXue et akl 
120111) . Stellar metalli cities, [Fe/Hl, for this sample are also 
available. Following IXue et all d201 ll) . we use the metal- 
licities obtained by the IWilhelm et al] meth odology in the 
SEGUE Stellar Parameter Pipeline (SSPP; see lLee et al.l2008l 
for details), which are likely to be the most reliable for stars 
with effective temperatures of BHB stars (on the order of 0.3 
dex). We then divide this sample on metallicity - 994 of 
our stars are relatively metal-rich ([Fe/H] > -2.0), while 871 
stars are very metal-poor ([Fe/H] < -2.0)0 

3. ANALYSIS METHOD 

3.1. Derivation of Rotational Velocity 

Let us denote by S an imaginary observer located at the Sun 
who is at rest with respect to the Galactic rest frame, and by O 
an observer located at the Sun who moves with the Sun. Here, 
we assume that the velocity of O with respect to S, Vq, and 
the three-dimensional (3D) position of the Sun with respect to 
the Galactic center are known. 

Now, suppose that the k-th star {k = 1 , • • • , n) is observed 
in the direction of x' k os by S and O. Then, the line-of-sight 
velocity with respect to S can be expressed as: 



Jos = v . „/os 



(1) 



where \ k is the velocity of this star with respect to S. Since 
the line-of-sight velocity of this star with respect to O is: 



los.kel / \ Ins 



we can calculate v l k os from (\q, v k ° s ' hel ,x' k os ). 
On the other hand, if we decompose as: 



Vfc = V,- k£r,k + Ve,k^9,k + V4>,ke<j>.k , 



(2) 



(3) 



we obtain 



where 



Jos . 



Q itk =x' k os -e Lk (i=r,9,(t > ). 



(5) 



Here, v,yt and e,^ (i = r,9,cjj) are the i-th velocity component 
of the A:-th star, and the basis unit vectors of the spherical co- 
ordinate system at the position of the k-th star, respectively. 
Given the assumptions: 

7 Our boundary metallicity ([Fe/H] = —2.0) lies between the peak metal- 
licities of the inne r-halo ([Fe/HJpeak ~ —1.6) and outer-halo ([Fe/HJpeak ~ 
—2.2) components (Carollo et al. 2007). 



• (Al) There is no correlation between the velocity and 
position of a star 

• (A2) The distributions of v r , vg, and v^ are symmetric 
around v r = Q,vg = 0, and = V rot , respectively 

then the data points {{Q^. kl v k m )} are likely to be distributed 
around a linear function of form: 



lt los 



(6) 



Thus, by performing a linear fit to the data in the Q$ - v los 
plane, we can obtain V rot - by which we denote the estimated 
value of V rot - by measuring the slope of th e best-fit linear 
functio n. This approach is similar to that of iFrenk & Whitel 
(fl980h . 

3.2. Derivation of 3D Velocity Dispersion 

In estimating the velocity di spersion, we use a modified ver- 
sion of the approach used by IWoolleyl (119781) . In principle, 
this method utilizes the direction dependence of the line-of- 
sight velocity dispersion. 

In addition to (Al) and (A2) above, let us further assume 
that the velocity distribution of the sample stars satisfies the 
following: 

• (A3) The velocity ellipsoid is aligned with a spherical 
coordinate system 

• (A4) The dispersions in the distributions of v r , vg, and 
v<£ around their centers are a>, og, and a^, respectively 

Then, it can be shown, for i = r,9,<p (see lMorrison et alll 19901 
for the case of i = d>), that 



E[var[v te e ; ]] 

" QtkQlk + rt 1 - E QlkQlk+vlH QtkQ 



n 



2 
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where we denote the expectation and variance of X by E[X] 
and var[X], respectively. 
Thus, by substituting the observed (var [v iM <2j] ) for 

E [var [v' os gi] ] > an d the already estimated V mt for V mt , we ob- 
tain: 
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Or 



(8) 



The solutions (of , a\ , a^) for this equation are the unbiased 
estimates for (a^,al,a^). Note that these estimated values 
are not guaranteed to take on positive values, so that we have 
to confirm the robustness of our method via an independent 
Monte Carlo simulation, as described below. 
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3.3. Testing the Reliability of our Method 

In order to estimate the error that accompanies our esti- 
mate of the three-dimensional velocity dispersion, we con- 
struct mock catalogs, and perform simulated observation of 
mock stars drawn from these catalogs. In this Monte Carlo 
simulation, each mock catalog is designed so that the distri- 
bution of r for the relatively metal-rich (or very metal-poor) 
mock stars resemble those of the observed relatively metal- 
rich (or very metal-poor) BHB stars, and that the velocity of 
the mock stars obey a given anisotropic Gaussian velocity dis- 
tribution (with or without net rotation). We vary the velocity 
anisotropy parameter, /3 = 1 — (p\ + a^) / (2of ), and produce 
1000 mock catalogs for each case. We find that the derived 
ag and a<f, is only reliable for r < 16 kpc and r < 18 kpc for 
the relatively metal-rich and very metal-poor sample, respec- 
tively, while estimates of o> are reliable at any r. 

4. APPLICATION OF OUR METHODOLOGY - DERIVATION OF 
ROTATION VELOCITIES, VELOCITY DISPERSIONS, AND 
VELOCITY ANISOTROPY PARAMETERS 

Throughout this paper, we assume that the Local Standard 
of Rest (LSR) is on a circular orbit with a rotation speed 
of 220 km s" 1 dKerr & Lvnden-Bell|[T986T) . It is worth not- 
ing that our assumed values for the Galactocentric distance 
of the Sun, Rq = 8.5 kpc, and the circular velocity of the 
LSR are both consistent with two recent inde p enden t de- 
terminations of these quantities by iGhez et al.l (120081) and 
iKoposov et ail ( f20loh . (Bow et alJ 42012ft have recently de- 
termined, on the basis of accurate line-of-sight velocities for 
stars in the APOGEE sub-survey of SDSS-III, that the circular 
velocity of the LSR is close to 220 km s" 1 . We also assume 
that the peculiar motion of the Sun with respect to the LSR 
jsJU Q ,V Q ,W Q ) = (10.0,5.3,7.2) km s" 1 dDehnen & Binnevl 
fl998h . 

FigureQ]shows the determination of mean rotational veloc- 
ity, y rot , for the relatively metal-rich (red) and very metal-poor 
(blue) BHB samples, as a function of Galactocentric distance, 
r. The shaded regions indicate the uncertainties in each result, 
estimated from a bootstrap approach (sampling with replace- 
ment). Inspection of this figure suggests that V lot at 13-23 kpc 
is a slightly decreasing function of r for the relatively metal- 
rich sample, while that for the very metal-poor sample is 
more or less flat. The relatively more metal-rich sample is 
in modest prograde rotation, with V, ot ~ 0-20 km s -1 , while 
the very metal-poor sample is in retrograde motion, with 
y rot <~ -20 to -50 km s" 1 , and lags that of the relatively metal- 
rich sample by ~ 20 -50 km s -1 , similar to prev ious observa- 
tions of stars in a much more local sample by ICarollo et al.l 
( 120071 120101) . t he outer-halo analys is based on distribution- 
function fitting (De ason et al.ll201 1[). and som e recently sim- 
ulated stellar haloes (e.g.. iTissera et al.l 12012ft . It is also in- 
triguing to see that the rotational shear between two samples 
seems to shrink at r ~ 20 kpc. We confirm that varying the 
LSR velocity only changes the absolute value of V rot , and that 
the very metal-poor sample always lags behind the relatively 
metal-rich sample, independent of the assumed LSR velocity. 

Figure|2]shows determinations of the radial velocity disper- 
sion, o>, and the two components of the tangential velocity 
dispersion, a§ and cr^, for the relatively metal-rich and very 
metal-poor samples, as a function of r. As seen in the fig- 
ure, a r is about 100 km s -1 , and exhibits a declining behavior 
over 12.5 < r < 20 kpc for both samples, supporting m ost 
previous studies (e.g. JXue et al . 2008; Brown et aT1l20 1 Oft . In 
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FIG. 1 . — The mean rotational velocity, V rot , for the relatively metal-rich 
BHB sample (red) and very metal-poor BHB samples (blue), as a function of 
Galactocentric distance r. Open circles show the result for our BHB samples. 
In this plot, we bin the relatively metal-rich (and very metal-poor) BHB stars 
in r by binning 500 stars (400 stars), sorted in r, and moving through the 
sample in steps of 20 stars. Each bin contains stars with a typical standard 
deviation in r of ~ 3 kpc, and the resultant V rot is presented at the median 
value of r. The associated shaded regions represent the uncertainties of our 
results, estimated from the bootstrap method, and denote the range covered 
by the 16% and 84% percentiles. 

addition, this figure suggests that <jg ~ holds for the very 
metal-poor sample. Noting that the gravitational potential is 
nearly spherical well above the disk plane, and that V rot of the 
very metal-poor stars is small compared with its velocity dis- 
persion, it follows that they obey a distribution function that 
depends mainly on the orbi tal energy and angular momentum 
dBinney & Trem aine 2008), and thus their spatial distribution 
is nearly spherical. The different behavior of the 3D veloc- 
ity dispersions for the two samples implies that both the func- 
tional form of the distribution function and spatial distribution 
for these samples are different. 

Figure [3] shows the velocity anisotropy parameter, j3 = 
1 — (<Tg + <r?)/(2of). This parameter quantifies the relative 
dominance of the radial and tangential velocity dispersions, 
and provides a simple diagnostic of the orbital properties of 
our sample stars. If j3 takes on values j3 < 0, we can in- 
fer that the halo stars are dominated by round orbits, while 
if < /? < 1, we can infer that the halo stars are dominated 
by radial orbits. When j3 is (nearly) zero, we can infer that 
the velocity distribution is (nearly) isotropic. Figure [3] sug- 
gests that relatively metal-rich BHB stars are dominated by 
radial orbits at 12 < r/kpc < 15 (/? = 0.3 ±0.3), while the 
very metal-poor BHB stars are dominated by circular orbits 
over 13 < r/kpc < 18 (/? = -0.9 ±0.7). A similar trend in 
the eccentricity distribution is also c onfirmed in the more lo- 
cal sample of halo stars reported by ICarollo et al.l (12010ft . in 
which very metal-poor halo stars possess a higher fraction of 
low-eccentricity orbits (see their Figure 5). 

5. DISCUSSION AND CONCLUSIONS 

The observed systematic differences in the orbital motions 
of the relatively metal-rich and very metal-poor BHB stars 
suggest that the orbital motion of the progenitor systems of the 
stellar halo depends on their metal abundances. Noting that 
dwarf galaxies with small er total stellar mas ses tend to have 
lower metal abundances dKirby et alJ 12008. and references 
therein), which can be understood as a result of a smaller 
rate of metal-enrichment events such as supernova explosions, 
our findings suggest that low-gas-fraction ('star-rich') sys- 
tems tend to move in radial orbits, while high-gas-fraction 
('gas-rich') systems tend to move in round orbits when they 
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FIG. 2. — Three-dimensional velocity dispersions o> (red), ag (blue), and 
cr^ (green), as a function of Galactocentric distance. Open circles show the 
result for our BHB samples. The same binning procedure as in Figure |T]is 
adopted. For the tangential velocity dispersions, some of the bins with large 
Galactocentric distance are excluded, because of a large systematic error indi- 
cated from Monte Carlo simulations. The associated shaded regions represent 
the uncertainties of our results, estimated from the bootstrap method. Results 
for the relatively metal-rich BHB sample and very metal-poor BHB sample 
are shown in panel (a) and (b), respectively. Filled squares at r = 8.5 kpc 
in panel (a) and (b) represent the Solar-neighborhood observations of halo 
stars with —2 < [Fe/H ] < —1.6 and [Fe/H] < -2, respectively, taken from 
IChiba & Yoshii 1 1993). 
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FIG. 3. — The velocity anisotropy parameter, /3 = 1 — (<r^ +cr^)/(2(T^), for 
the relatively metal-rich (red) and very metal-poor (blue) BHB samples, as 
a function of Galactocentric distance. Open circles show the result for our 
BHB samples, and correspond to those bins presented in Figure f2] The as- 
sociated shaded regions represent the uncertainties of our results, estimated 
from the bootstrap method. Filled squares at r = 8.5 kpc represent the Solar- 
neighborhood observations of h alo stars with —2 < [Fe /H] < —1.6 (red) and 
[Fe/H] < -2 (blue), taken from lChiba & Yoshiil fl99^) . The black solid and 
black dashed l i ne rep rese nt the f) profiles o f simulated stellar haloes from 
Diemand et al. 12005) and Sales et al. 12007), respectively. The dotted hori- 
zontal line at /? = is added to guide the eye. 



convert their gas to stars. 

In the hierarchical galaxy formation scenario, the Milky 
Way (and other large systems) attains it s mass as a result of 
the mergers of infalling smaller systems (IWhite & Reeslll978[ 
iBlumenthal et al.|[l984l) . Due to the deep gravitational poten- 
tial well of the M ilky Way, these in falling systems tend to 
have radial orbits (ISales et al.l 120071) . At this stage, we ex- 
pect that the orbital properties of infalling systems is indepen- 
dent of whether they are star-rich or gas-rich. However, the 
situation may be different when such radially infalling sys- 
tems pass near the Galactic center. In such a region, star-rich 
systems (e.g., massive dwarf galaxies with low gas fractions) 
are expected to be disrupted by tidal inter actions, which giv e 
rise to field halo stars with radial orbits (ISales et al.l 120071) . 
On the other hand, gas-rich systems are expected to interact 
with other gas-rich syst ems, and lose some orbital energy via 
dissipational processes (|Spmmer-Larsen & Christensenll989t 
lTheisin~996HSharma et al.l2012l) . If the angular momentum of 
such a gas-rich system is approximately conserved, the orbit 
is circularized, and the pericentric distance (distance of clos- 
est approach to the Galactic center) increases. Once its peri- 
centric distance become sufficiently large to avoid the central 
region of the Milky Way, further orbital change becomes less 
likely, since encounters with other gas-rich systems becomes 
less probable. Therefore, gas-rich systems tend to move in 
round orbits. When the Galactic disk forms, some gas-rich 
systems with similar orbital motions to disk gas may be ab- 
sorbed into the disk, due to their small relative velocities, 
while others remain moving on orbits across the halo, and 
eventually form field halo stars with round orbits. If the to- 
tal orbital angular momentum of gaseous systems is initially 
near zero, these halo stars would exhibit a net retrograde ro- 
tation with respect to disk stars. Some gas-rich systems in 
the Milky Way that have not yet formed many sta rs might be 
associated with the observ ed high-velocity clouds dBlitz et al. 
[T^IPulman et al-OoLl . 

In the Solar neighborhood, observations suggest that (3 ~ 
0.5 for halo stars, alm ost ind ependent of metal a bundance 
(IChiba & Yoshii|[l998T: see also lYoshii & Saiol[T979h . as over- 
plotted in Figure [3] Combined with our result, we see that 
the motion of the relatively metal-rich BHB stars move on 
radial orbits at 8.5 < r/kpc < 15. This behavior is consis- 
tent with recent cosmological simulations in which most halo 
stars originate from ac creted dwarf galaxies (iDiemand et alj 
l2005USales et alj2007h . as also overplotted in Figure|3] How- 
ever, our very metal-poor halo BHB stars suggest a transi- 
tion from radially-anisotropic to tangentially-anisotropic ve- 
locity distributions. The existence of very metal-poor outer- 
halo stars with round orbits, which is not confirmed in simu- 
lated stellar haloes, suggests that current simulations of disk 
galaxy formation may lack some important mechanisms, such 
as those proposed in the previous paragraph. We here note 
that many authors have pointed out that the simulated dwarf 
galaxies show an over-produ ction of stars, when compared 
with observed dwarf galaxies (Sa wala et al.ll201 lb . This phe- 
nomenon is often called the "overcooling problem", and this 
might result in underestimation of the gas-rich progenitors in 
simulated haloes, which in turn underestimate the numbers of 
metal-poor halo stars with round orbits. 

We find that the kinematics of outer-halo stars exhibit a 
marked dependence on stellar metal abundance, which pro- 
vides information about the physical properties of their pro- 
genitor systems (such as the gas fraction). Both of our main 
results, that the relatively metal-rich and very metal-poor stars 
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are dominated by radial and round orbits, respectively, and 
that the mean rotational velocity of very metal-poor halo stars 
lags that of relatively metal-rich halo stars, can be explained 
if very metal-poor stars originate from gas-rich systems and 
metal-rich stars from star-rich systems. Our findings cast a 
new light on the formation mechanism of the Milky Way and 



similar disk galaxies. 
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